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ABSTRACT 

This  report  presents  a  brief  summary  of  the  past  progress,  current  status 
and  development  trends  of  transonic  aerodynamics .  The  main  methods  to  com- 
pute steady  and  unsteady  transonic  flows  are  reviewed.  Also,  recent  advances 
in  transonic  buffet  prediction  and  transonic  airfoil  design  are  summarized 
and  attention  is  drawn  to  problems  requiring  further  intensive  research 
efforts . 


FOREWORD 

This  review  was  prepared  in  response  to  a  request  by  Mr.  Ray 
Siewert,  AIR-320,  Naval  Air  Systems  Command,  Washington,  D.  C. 

It  is  an  attempt  at  an  "overall"  summary  of  the  current  status  of 
transonic  aerodynamics  as  obtained  from  a  survey  of  the  published 
literature  and  from  discussions  with  selected  individuals.   Its  main 
purpose  is  to  provide  a  guide  to  current  work  and  development  trends  and 
to  draw  attention  to  problems  requiring  further  intensive  research  efforts 

The  author  gratefully  acknowledges  valuable  discussions  with  Mr.  Ray 
Siewert,  NASC,  Dr.  Zonars ,  Major  Butkewicz  and  Mr.  J.  Olsen,  AFFDL,  and 
Dr.  Yoshihara,  General  Dynamics. 

References  with  decimal  point  classification  refer  to  the  recent 
NASA  "Annotated  Bibliography  on  Transonic  Flow  Theory,"  NASA  TM  X-2363, 
September  1971.  All  other  references  are  listed  in  the  section 
"References." 
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TRANSONIC  AERODYNAMICS 

The  major  milestones  and  accomplishments  can  perhaps  briefly  be 
summarized  as  follows: 
19M+   Oswatitsch  (ref.  k.kSk) 

Formulation  of  transonic  small  perturbation  equation 
19^7   von  Karman  and  Oswatitsch  (ref.  J+.638  and  ref.  l) 

Formulation  of  transonic  similarity  rules  for  two-dimensional 

flow 
19^8   Lin,  Reissner  and  Tsien  (ref.  2) 

Classification  of  the  various  cases  and  criteria  for  linearization 
1950   Oswatitsch  (ref.  k.k6l   and  I+.U62) 

Development  of  "Integral  Equation  Method"  for  steady  two- 
dimensional  transonic  flow 
1950   Guderley  and  Yoshihara  (ref.  h.212) 

Development  of  "Hodograph  Method"  for  sonic  flow  past  rhombus 

profile 
1952   Oswatitsch  (ref.  k.h^k) 

Development  of  the  equivalence  rule 
1952   Whitcomb  (ref.  3) 

Development  of  the  area  rule 
1955   Oswatitsch  (ref.  ^57) 

Development  of  the  "Parabolic  Method"  for  sonic  flow  past  bodies 

of  revolution  and  airfoils 
1958   Spreiter  (ref.  U.568) 

Development  of  the  "Local  Linearization  Method"  for  sonic  flow 

past  airfoils  and  bodies  of  revolution 


i960   Sinnott  (ref.  2.3*0 

Development  of  a  "Semi -Empirical  Theory"  for  transonic  flow  past 
airfoils 

1961  Landahl  (ref.  k) 

Development  of  Linearized  Unsteady  Transonic  Flow  Theory 

1962  Guderley  (ref.  1.8) 

Publication  of  comprehensive  book  on  hodograph  solutions  for 

transonic  flows 
1962   Oswatitsch  (ref.  3.1) 

Organization  of  "Symposium  Transsonicum"  (first  international. 

conference  on  transonic  aerodynamics) 
1966   Whitcomb  (ref.  5) 

Development  of  supercritical  airfoil 

1966  Thomas  (ref.  6) 

Development  of  a  buffet  prediction  computer  program 

1967  Nieuwland  (ref.  U.U28) 

Development  of  hodograph  solution  for  shockless  transonic  flow 
past  quasi-elliptic  airfoils 

1968  Magnus,  Yoshihara,  MacKenzie,  Moretti,  Singleton  (ref.  3.2) 
Development  of  unsteady  finite  difference  method  for  two- 
dimensional  transonic  flow 

I968   AGARD  (ref.  3-2) 

International  specialists'  meeting  on  transonic  aerodynamics 
1970   Garabedian  &  Korn,  Murman  &  Cole,  Steger  &  Lomax,  Tai,  Norstrud 

(refs.  7,  U.l+23,  ^.581,  1+.591+,  h.hk5) 

Development  of  further   computer  programs   for  two-dimensional 

transonic  flow 


1972   NASA  (ref.  5) 

"Supercritical  Wing  Technology,"  Progress  Report  on  T-2C  and  F-8 
Flight  Evaluations 


BRIEF  CHARACTERIZATION  OF  THE  MAJOR  DEVELOPMENTS 

A.   Steady  Transonic  Flow  Theory 

During  the  period  from  World  War  II  to  the  Symposium  Transsonicum 
(I962)  two  main  approaches  were  developed  to  solve  the  nonlinear 
transonic  small  perturbation  equation,  i.e.  the  hodograph  method  and 
approximate  "direct"  methods. 

As  is  well  known,  for  two-dimensional  flow  the  hodograph  method 
leads  to  a  linear  equation  in  the  hodograph  plane.  However,  the 
satisfaction  of  the  flow  boundary  condition  becomes  quite  difficult  for 
general  airfoil  shapes.   Therefore,  solutions  could  be  developed  only 
for  special  configurations,  as  e.g.  for  the  rhombus  profile  (Guderley 
and  Yoshihara,  1950,  ref.  4.212).  These  results  as  well  as  a  thorough 
discussion  of  the  hodograph  method  are  given  in  Guderley' s  book  on  "The 
Theory  of  Transonic  Flow"  (1962). 

The  second  approach  (direct  methods)  was  pioneered  by  Oswatitsch 
(ref.  4.462)  who  recast  the  transonic  small  perturbation  equation  into 
an  equivalent  nonlinear  integral  equation  and  solved  it  approximately  by 
assuming  a  certain  decay  behavior  in  the  transverse  flow  direction. 
This  method  was  further  developed  by  Oswatitsch' s  student  Gullstrand 
(ref.  4.218)  and  by  Spreiter  at  NACA-Ames  (ref.  4.571).  An  alternate 
approximate  theory  was  also  introduced  by  Oswatitsch  (ref.  4.457) >  i.e. 
the  "parabolic  method."   In  this  approach  he  showed  that  for  half -bodies 
of  revolution  with  parabolic  meridian  profile  a  good  approximation  could 
be  obtained  by  linearizing  the  nonlinear  term  of  the  transonic  small 
perturbation  equation  and  thereby  reducing  it  to  a  parabolic  partial 
differential  equation.   Maeder  and  Thommen  (ref.  4.380)  extended  this 
method  to  study  transonic  flow  past  airfoils.   Spreiter  (ref.  4.568) 


showed  that  surprisingly  good  agreement  with  experiments  could  be 
achieved  by  applying  Oswatitsch's  linearization  only  "locally"  (the  local 
linearization  technique).  Another  extension  of  Oswatitsch's  parabolic 
method  was  proposed  by  Hosokawa  (nonlinear  correction  technique,  ref . 
U.253)*   In  this  method  a  correction  function  is  superimposed  upon  the 
linearized  Oswatitsch  solution  which  satisfies  a  nonlinear  ordinary 
differential  equation.   Its  solution  leads  to  the  prediction  of  transonic 
shocks.   Still  another  variation  of  Oswatitsch's  parabolic  method  is  due 
to  Cole  (ref.  U.106)  who  linearized  the  transonic  small  perturbation 
equation  in  such  a  way  that  the  mixed  elliptic -hyperbolic  character  of 
the  governing  equation  was  preserved. 

Brief  Summary  of  the  Work  Between  World  War  II  and  the  Symposium 
Trans sonicum  (1962):  Various  approximation  theories  were  developed  (the 
integral  equation  method,  the  parabolic  method,  the  local  linearization 
technique,  the  nonlinear  correction  method)  to  describe  two-dimensional 
and  axisymmetric  transonic  flow.  These  methods  are  valid  and  give  good 
agreement  with  available  experiments  only  over  a  limited  Mach  number 
range  and  for  special  categories  of  airfoils  and  bodies.   Only  the 
integral  equation  and  the  nonlinear  correction  technique  are  capable  of 
predicting  transonic  shocks.  Because  of  the  limitations  of  these  purely 
theoretical  methods  semi-empirical  approaches  must  be  relied  upon  for 
transonic  airfoil  design.   Sinnott's  method  becomes  the  generally  used 
tool  in  this  regard. 

During  the  last  decade  the  potential  of  the  high-speed  computer  for 
the  solution  of  transonic  flow  problems  became  more  and  more  apparent. 
By  1968  three  numerical  procedures  were  presented  at  the  AGARD 
symposium,  which  all  solved  the  two-dimensional  transonic  airfoil  case 


by  "brute  force,"  i.e.  the  steady  flow  was  obtained  as  the  asymptotic 
flow  for  large  times  in  an  unsteady  formulation,  where  the  unsteady 
Euler  equations  were  solved  by  a  proper  finite  difference  scheme  (Lax- 
Wendroff  scheme).  This  approach  has  the  advantage  that  the  mixed  elliptic- 
hyperbolic  character  of  the  steady-state  equations  is  replaced  by  a  purely 
hyperbolic  problem,  for  which  the  initial  value  problem  is  properly  set 
and  resulting  shocks  are  obtained  naturally  as  part  of  the  marching 
procedure.   This  approach  was  successfully  demonstrated  by  Magnus  and 
Yoshihara  (ref.  3.2),  MacKenzie  and  Moretti  (ref.  3.2)  and  Singleton 
(ref.  3.2). 

At  about  the  same  time  Nieuwland  of  the  National  Aerospace  Laboratory, 
Amsterdam,  was  able  to  obtain  shockless  transonic  potential  flow  solutions 
for  a  family  of  both  nonlifting  and  lifting  quasi-elliptical  airfoils 
using  the  hodograph  method.   This  work  not  only  supplied  a  reference 
solution  against  which  to  compare  the  above  described  numerical  approach 
but  also  clarified  the  long-standing  "transonic  controversy."  Furthermore, 
this  solution  gave  theoretical  support  to  the  experimental  studies  at  the 
National  Physical  Laboratory,  England,  by  Pearcey  and  collaborators. 
Their  work  was  based  on  the  idea  that  the  transonic  shock  strength  should 
be  reduced  by  a  special  design  of  the  airfoil  shape.  By  developing 
airfoils  with  a  "peaky"  pressure  distribution  Pearcey  indeed  succeeded  to 
virtually  eliminate  the  shock  wave  and  thus  to  reduce  the  transonic  drag 
rise. 

Nieuwland' s  work  in  turn  stimulated  further  use  of  the  hodograph 
equations.   In  combination  with  the  concept  of  "complex"  characteristics 
Garabedian  and  Korn  (ref.  7)  were  able  to  calculate  supercritical  airfoil 
sections  which  are  free  of  shocks  at  a  specified  speed  and  angle  of  attack. 


A  third  approach  was  developed  by  Murman  and  Cole  (ref .  U.U23)  who 
based  their  analysis  on  the  steady  transonic  small  perturbation  equation. 
Using  separate  difference  formulas  in  the  elliptic  and  hyperbolic  regions 
to  account  properly  for  the  local  domain  of  dependence  of  the  differential 
equation  the  transonic  potential  equation  was  solved  numerically. 

Garabedian  and  Korn  (ref.  9)  in  turn  developed  a  second  order 
accurate  version  of  the  Murman -Cole  method  and  also  incorporated  a 
conformal  mapping  procedure  (mapping  the  interior  of  the  unit  circle 
conformally  onto  the  exterior  of  the  profile)  thus  leading  to  a  quite 
desirable  distribution  of  mesh  points.  Work  along  similar  lines  was  also 
published  by  Jameson  (ref.  10). 

Another  refinement  of  the  Murman-Cole  relaxation  method  was  recently 
given  by  Steger  and  Lomax  (ref.  U.581)  by  using  the  full  nonlinear 
potential  equation  as  the  basic  equation  rather  than  the  small  disturbance 
equation. 

Also,  work  on  three  moi*e  methods  must  be  mentioned.  Tai  (ref.  h.^k) 
applied  the  method  of  integral  relations  to  the  steady  two-dimensional 
transonic  flow  problem.  Norstrud  (ref.  *+.U^5)  further  developed  the 
integral  equation  method  and  extended  it  to  compute  nonlinear  three- 
dimensional  transonic  flows  (ref.  U.W+).  Spreiter  (ref.  U.559)  applied 
his  method  of  local  linearization  to  the  computation  of  sonic  flow  past 
slender  bodies  and  wing-body  combinations.  Bailey  and  Steger  (ref.  11) 
and  Ballhaus  and  Bailey  (ref.  12)  extended  the  previously  described 
relaxation  technique  to  the  three-dimensional  transonic  flow  problem  and 
presented  results  and  comparisons  with  experiment  for  a  C-lUl  swept  wing 
model . 


Brief  Summary  of  the  Theoretical  Work  on  Steady  Transonic  Flows  During 
the  Past  Decade:   Spectacular  progress  has  been  made  in  our  ability  to 
compute  transonic  flows  with  the  introduction  of  large  capacity  high- 
speed computers.  The  following  major  developments  stand  out: 

1.  The  computation  of  shockless  two-dimensional  flows  by  hodograph 
methods  (Nieuwland,  Garabedian  and  Korn) . 

2.  The  computation  of  shocked  two-dimensional  flows  by  an  unsteady 
finite  difference  solution  of  the  Euler  equations  (Magnus  and  Yoshihara, 
MacKenzie  and  Moretti,  Singleton). 

3.  The  computation  of  two-  and  three-dimensional  transonic  flows  by 
relaxation  methods  (Murman,  Cole,  Steger,  Lomax,  Garabedian,  Korn, 
Jameson,  Bailey,  Ballhaus ,  Krupp) . 

h.     The  computation  of  three-dimensional  transonic  flows  by  the  integral 
equation  method  and  the  method  of  local  linearization  (Norstrud, 
Spreiter) . 

5.  Recent  comprehensive  reviews  of  these  developments  have  been  given  by 
Yoshihara  in  refs.  19  and  36. 
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B.   Unsteady  Transonic  Flow 

In  19^8  Lin,  Reissner  and  Tsien  (ref.  2)  showed  that  linearization 
of  the  governing  equations  is  possible  whenever  the  flow  is  "sufficiently" 
unsteady.  Hence,  the  computation  of  transonic  flow  past  oscillating  wings -- 
as  required  for  flutter  and  dynamic  stability  calculations --can  be  based  on 
the  linearized  unsteady  potential  equation  provided  the  condition 

k  »  (1  -  ML) 

is  satisfied  (k  =  reduced  frequency,  1VL  =  local  Mach  number). 

This  linearized  equation  was  studied  in  detail  by  Landahl  (ref.  k) 
who  obtained  solutions  for  oscillating  airfoils,  low  aspect  ratio  wings, 
low  aspect  ratio  wing-body  combinations,  semi-infinite  rectangular  wings, 
rectangular  wings  of  arbitrary  aspect  ratio  and  wings  of  convex  and  concave 
polygonal  planforms.  A  comprehensive  review  of  this  work  was  given  by 
Landahl  in  his  book,  "Unsteady  Transonic  Flow,"  ref.  k.     The  need  for  a 
generally  applicable  "lifting  surface"  method  led  Watkins  and  collaborators 
at  NACA-Langley  to  the  development  of  the  "kernel  function  method"  for 
oscillating  finite  wings  of  arbitrary  planform  in  subsonic  flow  (ref.  8). 
Its  extension  to  sonic  flow  was  given  by  Runyan  and  Woolston  in  1957 
(ref.  13).  Stark  of  Saab  Airplane  Company  (ref.  lU)  developed  a  similar 
lifting  surface  theory  for  sonic  flow  by  also  using  a  linear  approximation 
to  the  lift  distribution  but  determining  the  weight  coefficients  by 
satisfying  the  tangency  condition  in  a  least  square  sense.  Rodemich  and 
Andrew  (ref.  15 )  solved  this  same  problem  by  using  the  transonic  counter- 
part of  an  earlier  supersonic  numerical  lifting  surface  theory,  i.e.,  the 
"transonic  box  method."  In  a  further  study  Stenton  and  Andrew  treated  the 
transonic  aerodynamics  for  planar  wings  with  trailing  edge  control  surfaces, 


In  all  these  studies  an  infinitely  thin  lifting  surface  oscillating 
at  "sufficiently"  high  reduced  frequency  had  to  "be  assumed. 
Unfortunately,  in  many  practical  stability  and  flutter  problems  the 
reduced  frequencies  are  quite  low,  thus  invalidating  the  basic 
assumptions  of  linearized  unsteady  transonic  flow  theory.  Hence,  the 
aerodynamic  forces  on  slowly  oscillating  transonic  wings  must  be  found 
by  superimposing  small  unsteady  disturbances  upon  the  mixed  subsonic, 
supersonic  nonuniform  steady  transonic  flow  field.  Thus,  the  basic 
nonlinear  equation  must  be  used  which  now  represents  a  form  of  local 
wave  equation  wherein  the  local  convective  velocity  and  the  local  speed 
of  sound  must  be  determined  from  the  steady  flow  calculation.  The  use 
of  velocity  potentials  for  nonuniform  transonic  flow  was  first  suggested 
by  Landahl  and  an  approach  was  outlined  by  Andrew  and  Stenton  (ref.  U.ll). 
It  is  based  upon  the  application  of  Fermat's  principle,  i.e.,  the  methods 
of  geometric  acoustics  are  used  to  determine  acoustic  ray  paths  and 
transmission  times.  An  approximation  theory  for  nonlinear  unsteady 
transonic  flow  was  suggested  by  Teipel  (ref.  3.1)  and  Hosokawa  (ref.  3-1) > 
who  generalized  the  previously  discussed  "parabolic  method"  and  the 
"nonlinear  correction  theory"  to  sonic  flow  past  oscillating  airfoils. 
An  application  of  these  methods  to  transonic  flow  past  slowly  oscillating 
bodies  of  revolution  was  given  by  Liu,  Platzer  and  Ruo  (ref.  16). 

Summary :   The  linearized  unsteady  transonic  flow  theory  is  reasonably 
well  developed  for  planar  wings  of  arbitrary  planform.  However,  this 
theory  holds  only  for  infinitely  thin  wings  oscillating  at  sufficiently 
high  frequency.   Practical  stability  and  flutter  problems  generally  occur 
at  low  reduced  frequencies.   In  addition,  dynamic  pressure  for  flutter 
often  exhibits  a  significant  dip  near  M  =  1,  so  that  the  greatest  need 
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for  reliable  flutter  predictions  arises  at  transonic  flight  speeds.  No 
method  exists  at  the  present  time  to  satisfy  this  need. 
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C.  Transonic  Buffeting  Prediction 

The  only  method  attempting  to  predict  the  onset  of  transonic 
buffeting  developed  to  date  is  that  published  by  F.  Thomas  (ref.  6).   It 
is  restricted  to  two-dimensional  airfoils  and  requires 

1.  the  calculation  of  the  inviscid  supercritical  pressure  distribution 
including  the  position  of  the  shock. 

2.  'a  boundary  layer  calculation  to  determine  the  separation  point. 

3.  the  incorporation  of  a  reliable  buffet-onset  criterion. 

Thomas  defines  buffet  onset  to  be  reached  when  the  calculated  boundary- 
layer  separation  point  has  moved  forward  on  the  airfoil  so  that  it  . 
coincides  with  the  computed  shock  position.   Other  criteria  were  proposed 
by  Outman  and  Lambert  (ref.  17),  Gadd  (ref.  18),  Sinnott -Osborne  (ref.  18) 

A  detailed  review  of  the  available  inviscid  and  viscous  analysis 
methods,  boundary -layer  transition  and  separation  and  shock  boundary- 
layer  interaction  analyses  was  recently  given  by  Gentry  and  Oliver 
(ref.  20),  who  also  described,  tested  and  documented  the  Thomas  buffeting 
computer  program. 

This  program  uses  the  Sinnot-Osborne  semi -empirical  method  for  the 
prediction  of  the  transonic  pressure  distribution  and  the  shock  location. 
The  boundary-layer  calculation  is  based  on  an  integral  method  by  Walz 
(ref.  21).  Although  such  a  boundary-layer  calculation  cannot  account  for 
the  difficult  shock-boundary-layer  interaction  mechanism  good  agreement 
with  the  few  available  test  data  was  obtained. 

Summary:   The  availability  of  large  capacity  high-speed  computers  has 
made  it  possible  to  attempt  the  theoretical  prediction  of  transonic 
buffet  onset.   One  computer  program — originally  developed  by  F.  Thomas 
in  Germany—has  been  tested  and  checked  at  McDonnel  Douglas  and  is  now 
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available  for  general  use.  Also,  a  comprehensive  survey  of  the  transonic 
"buffet  problem  has  recently  been  completed  (ref.  20) .  The  present  computer 
program  is  valid  only  for  the  two-dimensional  case  and  contains  components 
of  doubtful  validity.  Buffet  onset  prediction  for  three-dimensional  wings 
is  still  unsolved. 
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D.  Transonic  Wing  Design 

The  major  achievements  to  date  have  been  fairly  well  documented 
(refs.  5,  2.23).  Therefore,  only  a  brief  review  of  the  major  ideas  and 
results  is  presented. 

The  transonic  design  objectives  are  well  known,  i.e., 

1.  for  the  cruise  condition  the  primary  goal  is  to  increase  the  drag 
divergence  Mach  number  to  as  large  a  value  as  possible  while 
maintaining  a  prescribed  lift  coefficient  and  buffet  margin. 

2.  for  the  maneuver  condition  the  primary  goal  is  to  delay  buffet  onset 
to  as  high  a  value  of  lift  coefficient  as  possible  for  a  given  Mach 
number. 

These  two  transonic  phenomena,  namely  drag  divergence  and  buffet, 
are  caused  by  the  formation  of  shock  waves  on  the  airfoil  upper  surface, 
i.e.,  the  drag  divergence  is  produced  by  the  shock  losses  (entropy 
increase)  and  the  buffeting  is  generated  by  shock  induced  boundary  layer 
and  leading-edge  separation.  Hence,  the  transonic  design  goal  consists 
of  the  development  of  wing  profiles  which  produce  only  a  weak  shock  while 
maintaining  a  prescribed  lift  coefficient  at  as  high  a  Mach  number  as 
possible. 

The  following  features  will  clearly  help  to  achieve  this  goal: 

1.  Lift  should  be  produced  as  far  as  possible  by  overpressures  on  the 
lower  surface  without  adversely  affecting  the  flow  over  the  upper 
surface . 

2.  Lift  from  the  upper  surface  should  be  obtained  by  extending  the 
chordwise  extent  of  the  underpressures  and  increasing  its  level  with- 
out increasing  the  shock  strength. 
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3.  Drag  will  be  minimized  if  underpressures  on  the  upper  surface  are 
acting  on  upstream  facing  surface  elements,  overpressures  on  the 
lower  surface  on  downstream  facing  surface  elements. 
These  considerations  have  helped  to  evolve  the  following  essential 

ideas  in  transonic  airfoil  design: 

1.  Pearcey's  "peaky"  airfoil  (ref.  2.31).  The  airfoil  is  designed  with 
a  large  leading  edge  radius  thereby  producing  expansion  waves  at  the 
nose  which  are  reflected  from  the  sonic  line  as  compression  waves 
thus  decreasing  the  shock  strength. 

2.  Whitcomb's  "aft  cambered"  airfoil  (ref.  5).  Aft  camber  produces  lift 
by  increasing  the  lower  surface  overpressures  and  by  extending  the 
upper  surface  underpressures  without  strengthening  the  upper  surface 
shock  wave. 

3.  Yoshihara's  "humped  upper  surface"  airfoil  (ref.  22).  The  Mach 
number  ahead  of  the  shock  is  lowered  by  "humping"  the  upper  surface 
(Yoshihara  introduces  the  concept  of  the  "quasi -limiting  Mach  wave" 
to  demonstrate  the  effectiveness  of  this  concept). 

k.      Jet -flapped  airfoil.  This  concept  was  first  proposed  by  Poisson- 
Quinton  (ref.  23).  The  effect  of  the  jet  flap  is  to  increase  the 
overpressures  over  the  entire  lower  surface  and  to  move  the  shock 
wave  toward  the  trailing-edge. 
Progress  in  the  development  of  successful  aft -cambered  airfoils  has 
been  spectacular.  The  wind  tunnel  and  flight  testing  (T-2C  and  F-8 
airplane)  have  demonstrated  the  superiority  of  the  "supercritical  airfoil." 
Its  performance  characteristics  are  fully  documented  in  ref.  5.  A  joint 
NASA/USAF  program  is  now  in  progress  which  will  generate  additional 
information  on  the  supercritical  wing  characteristics  of  an  F-11A 
variable  sweep  aircraft.  The  specific  objectives  of  this  TACT  program 

15 


(Transonic  Aircraft  Technology)  are  to  demonstrate  transonic 
maneuverability  by  evaluating  buffet  onset  and  intensity,  determining 
overall  performance  and  handling  qualities  and  measuring  local  pressure 
distributions  and  wake  drag. 

The  state  of  the  art  of  high  lift  devices  was  recently  reviewed  by 
Fairchild  Hiller  Corporation  (ref.  2k).     Although  some  of  the  statements 
made  in  this  report  concerning  the  present  status  of  theoretical 
techniques  at  transonic  speeds  have  been  obsoleted  by  the  rapid  develop- 
ment of  the  past  three  years,  a  theoretical  prediction  of  the  aerodynamic 
characteristics,  buffet  onset  and  established  buffet  on  airfoils  with 
high  lift  devices  is  still  impossible.  Furthermore,  sufficiently 
detailed  experimental  information  is  also  unavailable  to  establish  the 
buffet  onset  mechanism  of  leading-edge  devices.  A  similar  situation  is 
found  for  trailing-edge  devices.  However,  the  assumption  of  a  bubble- 
type  mechanism  is  probably  quite  correct.   Some  information  has  been 
obtained  since  the  publication  of  the  Fairchild-Hiller  report  on  jet 
flaps  indicating  quite  encouraging  buffet  improvements  (refs.  22,  25,  26, 
27). 


16 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Steady  Transonic  Flow  Theory 

#    Three  new  methods  have  been  developed  over  the  past  few  years 

which  make  it  possible  to  obtain  "exact"  solutions  for  inviscid 
steady  transonic  flow  over  two-dimensional  airfoils,  i.e., 

I.  the  hodograph  approach  for  shockless  transonic  flow 
(Nieuwland&  Boerstoel,  Garabedian  &  Korn). 

II.  the  unsteady  finite  difference  approach  (Magnus  &  Yoshihara, 
Moretti,  Singleton). 

III.  the  relaxation  method  (Murman  &  Cole,  Lomax  &  Steger, 
Bailey,  Ballhaus ,  Jameson). 

£    The  last  two  methods  are  the  most  useful  methods  for  general 
design  purposes. 

%  Comparing  these  two  methods  the  following  features  stand  out: 

The  unsteady  finite  difference  approach  requires  much  more 
computer  time.  The  unsteady  finite  difference  method  is  mathe- 
matically well  founded  (using  Lax's  concept  of  weak  solutions  of 
hyperbolic  initial  value  problems)  resulting  in  the  correct 
capture  of  the  shock. 

The  relaxation  method,  in  contrast,  does  not  properly 
capture  the  shock  and  further  work  is  necessary  to  remedy  this 
deficiency.  Potentially,  however,  the  relaxation  procedure 
(especially  when  applied  to  the  exact  potential  equation  rather 
than  to  the  small  disturbance  equation)  is  the  most  promising 
method  because  of  the  significantly  lower  computer  time  required 
and  because  of  its  ready  extendability  to  three-dimensional  flows 
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Three  methods  are  presently  under  development  to  compute  three- 
dimensional  transonic  flows,  i.e., 

I.  the  local  linearization  method  (Spreiter). 

II.  the  integral  equation  method  (Norstrud) . 

III.  the  relaxation  method  (Bailey,  Ballhaus,  Steger). 

The  local  linearization  method  is  restricted  to  flows  close 
to  M  =  1  and  cannot  predict  the  occurrence  of  shocks.  The 
relaxation  method  appears  as  the  more  promising  method  when 
compared  with  the  integral  equation  method  since  no  further 
approximations  need  to  be  introduced.  However,  the  previously 
mentioned  difficulty  of  correctly  capturing  the  shock  still  needs 
to  he  overcome. 

The  prediction  of  viscous  effects  still  remains  an  extremely 
difficult  problem.  The  development  of  simplified  transonic  shock- 
boundary-layer  models  remains  probably  the  only  realistic 
alternative  in  the  absence  of  rational  describing  equations. 
Yoshihara's  proposal  of  a  "viscous  ramp"  coupled  with  an 
experimental  determination  of  the  entrainment  parameters  as  well 
as  other  previously  proposed  models  (see  refs.  20  and  28  for 
recent  comprehensive  surveys)  should  form  reasonable  starting 
points  for  further  progress. 
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B.   Unsteady  Transonic  Flow  Theory 

#  Linearized  unsteady  transonic  flow  theory  is  fairly  well 
developed.  However,  this  theory  is  of  little  use  because  of  its 
high-frequency  limitation. 

%        No  theory  exists  at  the  present  time  to  cover  the  lower  frequency 
range  needed  for  practical  flutter  and  dynamic  stability 
calculations. 

#  NASA-Langley  is  presently  supporting  efforts  to  develop  such  a 
theory  using  relaxation  techniques  (Boeing),  the  local 
linearization  technique  (Nielsen  Engineering),  the  general 
aerodynamic  lifting  surface  element  approach  (General  Dynamics) 
and  the  method  of  geometrical  acoustics  (Lockheed) . 

%        All  these  methods  are  based  on  inviscid  flow  theory.  Unsteady 
shock-boundary-layer  interactions  therefore  are  excluded. 

#  Very  little  experimental  information  is  available.  A  few 
detailed  oscillatory  pressure  measurements  on  two-dimensional 
airfoils  have  been  obtained  in  the  Netherlands  and  Germany. 
Pressure  measurements  on  an  oscillating  delta  wing  have  recently 
been  completed  at  Cornell  Aero  Lab.  Much  more  data  of  this  type 
is  required  to  evaluate  the  theoretical  approaches. 
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C.  Transonic  Buffeting 

#  A  method  has  been  developed  (Thomas)  to  predict  the  onset  of 
transonic  buffeting  on  two-dimensional  airfoils. 

#  Further  refinements  and  detailed  comparisons  with  well- 
controlled  experiments  are  necessary. 

#  The  three-dimensional  problem  still  defies  rational  analysis. 
However,  the  recent  successful  development  of  three-dimensional 
inviscid  methods  in  conjunction  with  the  Thomas  program  could 
open  up  some  possibilities  in  the  near  future. 

#  A  thorough  review  of  this  problem  area  has  been  recently 
completed  under  ONR-sponsorship. 
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D.  Transonic  Wing  Design 

#  Spectacular  progress  has  been  made  in  the  development  of  two- 
dimensional  airfoils  suitable  for  operation  close  to  Mach  One 
(supercritical  airfoil). 

#  In  contrast  to  the  1968  status  assessment  (see  AGARD  Report  No. 
17,  ref .  2.20)  which  deplored  the  low  level  of  work  on  three- 
dimensional  problems,  both  theoretical  and  experimental  efforts 
(wind  tunnel  and  flight  test)  have  recently  been  intensified  and 
should  yield  valuable  design  information  (F-111A  variable  sweep 
flight  tests,  NASA-Ames  three-dimensional  relaxation  procedures, 
wind  tunnel  tests  of  wing -body  combinations,  ref.  29,  wind 
tunnel  tests  of  three-dimensional  jet-flapped  wing,  ref.  30). 

#  Little  improvement  has  been  made  in  the  development  of  reliable 
design  procedures  for  high-lift  devices  since  the  publication  of 
the  1969  comprehensive  review  (ref.  2^). 

9    However,  recent  studies  of  the  jet-flap  effect  shows  promising 
new  possibilities  to  improve  the  buffeting  characteristics. 
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As  a  result  of  the  above  outlined  recent  results  and  trends  the 
following  recommendations  are  suggested: 

#  Further  work  is  needed  on  the  development  of  faster  and  simpler, 
yet  sufficiently  accurate  computation  methods  for  steady, 
inviscid,  planar  transonic  flows.  Work  along  these  lines  is 
presently  being  carried  out  by  Magnus  and  Yoshihara  for  NASA- 
Ames,  ref.  31.  This  work  appears  to  be  the  most  advanced  and 
therefore  should  be  carefully  assessed  when  considering  research 
proposals  in  this  area. 

#  Vigorous  efforts  should  be  undertaken  to  incorporate  viscous 
effects  in  the  computation  of  steady,  planar  transonic  flows. 

%        A  systematic  experimental  program  should  be  formulated  and 

supported  to  provide  a  reliable  data  basa  for  comparison  with 
theory  (see  also  the  suggestions  in  ref.  20). 

9    Vigorous  support  should  be  given  to  work  on  three-dimensional 
flows,  wing -body  interference,  engine  nacelle  installation, 
viscosity  effects,  and  the  evolution  of  practical  design  criteria, 

%        Particular  emphasis  should  be  given  to  the  study  of  high-lift 

devices.  The  jet-flap  appears  to  offer  excellent  possibilities 
to  improve  the  maneuvering  characteristics.  Further  testing  of 
the  maneuvering  jet  flap  concept  is  required. 

#  The  development  of  transonic  dynamic  stability  and  flutter 
prediction  methods  should  be  vigorously  supported. 

#  The  measurement  of  detailed  transonic  oscillating  pressure 
distributions  should  be  initiated  and  supported.  The  NLR- 
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technique  (ref .  32)  which  requires  only  one  pressure  transducer 
might  offer  considerable  advantages.  (See  also  ref  s .  33  and  3^- 
for  experience  with  this  technique  in  the  United  States . ) 

Support  should  he  given  to  the  development  of  promising  new 
measuring  techniques .   Interf erometric  holography  seems  to  offer 
new  possibilities  as  exemplified  by  the  recent  work  of  d.  J. 
Collins,  ref.  35  and  should  be  further  used  to  provide  detailed 
information. 
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3   ABSTRACT 


This  report  presents  a  brief  summary  of  the  past  progress,  current  status  and 
development  trends  of  transonic  aerodynamics.  The  main  methods  to  compute  steady 
and  unsteady  transonic  flows  are  reviewed.  Also,  recent  advances  in  transonic 
buffet  prediction  and  transonic  airfoil  design  are  summarized  and  attention  is  drawn 
to  problems  requiring  further  intensive  research  efforts. 
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